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Abstract 



^ , In supersymmetric theories, the decays of the neutral CP-even and CP-odd as 

H I well as the charged Higgs bosons into scalar quarks, in particular into top and 

bottom squarks, can be dominant if they are kinematically allowed. We calculate 
the QCD corrections to these decay modes in the minimal supersymmetric extension 
of the Standard Model, including all quark mass terms and squark mixing. These 
corrections turn out to be rather large, altering the decay widths by an amount 
which can be larger than 50%. The corrections can be either positive or negative, 
and depend strongly on the mass of the gluino. We also discuss the QCD corrections 
to the decays of heavy scalar quarks into light scalar quarks and Higgs bosons. 



1. Introduction 

Supersymmetric theories [|l| are widely considered as the most attractive extensions of 
the Standard Model. In a grand unified framework, they protect scalar Higgs bosons 
from acquiring very large masses and provide the opportunity to generate the electroweak 
symmetry breaking radiatively. In the minimal version of these theories, the Minimal 
Supersymmetric Standard Model (MSSM), the Higgs sector is extended to contain two 
doublet fields, leading to five physical states 0: two CP-even neutral Higgs bosons, h 
and H, a CP-odd Higgs boson A and two charged Higgs particles H^. While the lightest 
CP-even Higgs boson h mass is predicted to be less than 130 GeV @, the H, A and H^ 
states are expected to have masses of order of the electroweak symmetry breaking scale. 

If all genuine SUSY particles are very heavy, the neutral and charged Higgs bosons 
will decay into standard fermions and gauge bosons, as well as into cascades involving 
the lighter Higgs bosons; in most cases the decays into heavy b and t quark pairs are 
dominant. These decay modes have been extensively discussed in the literature and 
for a recent summary we refer the reader to Ref. ^. However, it could well be that 
charginos, neutralinos and sfermions are light enough for the Higgs decays into these 
particles to be kinematically allowed. For instance, in grand unified models with proper 
radiative electroweak symmetry breaking, the H, A and H^ bosons are rather heavy 
and approximately mass degenerate, Mh — Ma — Mh± ~ a few hundred GeV, while 
charginos, neutralinos and eventually top and bottom squarks have masses of (9(100 
GeV) 1^. In this case, the decay pattern of the heavy Higgs particles H,A and H^ will 
be drastically different [^ |31- In particular, because of the large Yukawa couplings and 
large squark mixing, the decays of the heavy Higgs bosons into top and bottom squarks 
can be competitive with the standard decay channels and can even be the dominant ones 
in some areas of the MSSM parameter space 0. 

It is well known that the standard hadronic decay modes of Higgs particles are signifi- 
cantly affected by QCD radiative corrections p|. In order to have full control on the Higgs 
decay widths into squark pairs - and to make a reliable comparison with the standard 
decays - QCD corrections have also to be included. Recently, the QCD corrections to the 
decay of the charged Higgs boson into tb pairs have been calculated and found to be 
quite substantial. 

In the present paper, we derive the QCD corrections to the decays of the neutral 
CP-even and CP-odd Higgs bosons into stop and sbottom pairs. These corrections are 
found to be rather large, altering the Born decay widths by an amount which can exceed 
50%. These effects must therefore be taken into account. We also rederive the QCD 
corrections to the decays of the charged Higgs boson into tb pairs using a different [and 
more convenient] renormalization scheme compared to the one adopted in Ref. 0. We 
find that the corrections can be extremely large and can strongly suppress the decay 
width compared to its tree-level value. Finally, we also discuss for completeness the QCD 
corrections to the reverse decay of heavy squarks into their lighter partners plus light 
Higgs bosons; these decays have been discussed at tree-level in Ref. [l^ . 



The paper is organized as follows: to set the notation and to introduce the various 
parameters which enter the analysis, we summarize in the next section the Higgs decay 
widths into squark pairs in Born approximation. In section 3, we present the analytic 
expressions of the QCD corrections to the decay widths. Section 4 contains the discussion 
of the numerical results. In section 5, we discuss the QCD corrections to the decays of 
heavy squarks into Higgs bosons and we give our conclusions in section 6. 

2. Born Approximation 

The amplitudes for the decay widths of the MSSM neutral heavy CP-even, CP-odd and 
charged Higgs bosons[], that we will collectively denote by $, into the scalar partners of 
first and second generation quarks 

$ ^ g,g; , <^ = H,A,H^ (1) 

depend on three parameters if the quark mass is neglected: the mass of the decaying Higgs 
boson, the squark masses [more precisely the left- and right-handed soft SUSY breaking 
scalar masses m^^ and rrigj^, which in general are taken to be equal] and the ratio of the 
vacuum expectation values of the two MSSM Higgs fields tan/3. A mixing angle a in the 
CP-even Higgs sector also enters the amplitudes, but in the MSSM it can be expressed 
in terms of tan/5 and the Higgs boson mass M$. 

In the case of the third generation scalar quarks, and in particular top squarks, the 
mixing between left- and right-handed squarks, which is proportional to the mass of the 



partner quark, must be included |T^. In terms of the scalar masses rriq^^mq^^ the Higgs- 
higgsino mass parameter n and the soft SUSY-breaking trilinear coupling Aq, the squark 
mass matrices read 



m] + MIl mq Mlr 
'^ \ mqMLR ml + Ml^ 



m ' '''^ 

MIl = ml + {II -eqsl.) cos 2P Ml 

MIr = ml + eqs''^ cos 2/3 Ml 

Mlr = Aq + fi{t^n/3)-^'l . (2) 

/f and Cq are the weak isospin and electric charge of the squark q, and s^ = 1 — c^ = 
sin^^^. They are diagonalized by the 2x2 rotation matrices 

R^ = { ' ' , cg„ = cos On and sg„ = sin 9a (3) 



-'^In view of the experimental bounds on the squark masses |0|, the hghtest CP-even Higgs boson 
in the MSSM cannot decay into squark pairs since its maximal mass value is smaller than Af/,, < 130 
GeV P]. These decays will therefore not be considered here. However, the analytical expressions can be 
straightforwardly obtained from the decays of the heavy CP-even Higgs boson, after the proper changes 
in the Higgs boson couplings. 



which turn the mass eigenstates, gi and g2, into the current eigenstates ql and qR, the 
mixing angle and squark masses are then given by 



tan6'„ 



2 m„ Mlr 



Mh - Ml^ - J{Mh - Ml^f + 4 ml Ml^ 



(4) 



and 



"^1,2 = "^? + 2 



1 r 



ML + Ml^ T J {Ml, - Ml^f + WMl^ 



(5) 



In the Born approximation, the decay widths of the MSSM heavy CP-even, CP-odd 
and charged Higgs bosons into squark pairs, eq. ([|), can be written as [|13 






\'/'{Mlml,ml) iG^,,y 



(6) 



where \ = x"^ + y"^ + z^ — 2 {xy + xz + yz) is the two-body phase space function. The 
couphngs of the Higgs bosons to squarks G^ij read 



1 ^ - T 



(7) 



with the matrices Cc^qqi summarizing the couphngs of the Higgs bosons to the squark 
current eigenstates; for the if, A and H^ particles, they are given by 

_ / (2/| - 2eqsl,) Ml cos(/5 + a) + 2m2r? mq{Aqrl + ^irl) \ 

""'' \ mq{Aqr\ + ^irl) 2eqs'^MlcosiP + a) + 2my, ) ^^^ 



C 



-m„ 



Aqq 



m„ 



fi — Ag(tanP) 



-2/? 



/i — AgitanP) 




-2/| 



(9) 



C 



H±ib 



r- ( m^ tan/3 + m^/tan/5 — M^ sin 2/3 m^ (A^tan/? — /z) 
I vfit [At/ tan (3 — jj) 2mt vn^^j sin 2/3 



with the coefficients r\ 2 as 



(10) 



sma 
sin/3 



cos a 



sin/3 



cos a 



cos/3 



sma 



cos/? 



'Ill 



In principle, for the Higgs boson masses and the mixing angle a, one has also to include 
the large radiative corrections |^ which grow as the fourth power of mt and induce an 
additional dependence on jj, and Aq at the subleading level. However for rather heavy 
Higgs bosons, Mh ~ Ma ~ Mh± ^ M^, these corrections will not affect the decay 



amplitudes: the decaying Higgs boson mass can be used as input parameter and the 
angle a will reduce anyway to the value a — >• /3 — 7r/2. In this case, the couplings of the 
CP-even Higgs boson will simplify to 

^ / (2J|-2e,s2^.)M|sin2/?T2m2(tan/3)=Fi Tm,[A,(tan/?)Ti - ^] \ 

""^ \ Tmg[Ag(tan/?)Ti-/i] 2egS^M|sin2/5 ^ 2m^(tan/3)Ti j 

where the "-" is for up-type squarks and the "+" for down-type squarks. When the 
quark masses and the squark mixing angles are set to zero, as is the case for first and 
second generation squarks, the pseudoscalar Higgs boson cannot decay at tree-level into 
squark pairs since the AqiQi coupling is zero by virtue of CP-invariance and the Aqiq2 
coupling is proportional to mq. 

The decay widths of the heavy neutral CP-even and the charged Higgs bosons into 
first and second generation squarks are proportional to 

r° ~ GpM^ sin^ 2/3/M.j. 

in the asymptotic regime M$ ^ nig. These decays are suppressed by the heavy Higgs mass 
and therefore cannot compete with the dominant decay modes into top and/or bottom 
quarks [and to charginos and neutralinos] for which the decay widths grow as M$. 

In contrast, for the decays involving top squarks, the partial widths up to mixing angle 
factors are proportional to 

r° ~ Gprnf/iM^ tan^ /3) or/and Gprn^^ifi + At/ tan (3f/M^ . 

For small tan jS values and not too heavy Higgs bosons, or for intermediate values of tan (5 
and for fi and At values of the order of ~ 1 TeV, the partial decay widths into top squarks 
can be very large and can compete with, and even dominate over, the decay channels into 
top quarks [and into charginos/neutralinos] . Furthermore, decays into sbottom quarks can 
also be important for large values of tan/3 and A^. We therefore focus in this paper on the 
Higgs boson decays into third generation squarks; however, since we give complete analytic 
expressions, the decay widths into first and second generation can straightforwardly be 
obtained by setting the quark masses and the squark mixing angles to zero. 

3. QCD Corrections 

The QCD corrections to the Higgs decays into scalar quarks, eq. (|l]), consist of virtual 
corrections, Figs.la-d, and real corrections with an additional gluon emitted off the final 
squark states. Fig. le. The 0{as) virtual contributions can be split into contributions 
with gluon [la] and gluino [lb] exchange in the ^qq' vertex, the mixing contribution 
due to the off-diagonal self energies of the outgoing squarks [Id] and the quartic squark 
interaction contribution [Ic]. These contributions have to be supplemented by the wave- 
function counterterms, and by a counterterm renormalizing the $gjg' interaction. 



The partial widths of the decays of the MSSM Higgs bosons into scalar quarks, eq. (|l|), 
at 0{aas) can be written as 



r^($ 



Q^Qi) 



Gf 



a. 



2^/2txMI Ti 



A^/2(M|,m? ml) G$,, A 



(J)jj i-^<^ij 



with 



A 



<j>ij 



A y I ACT , A fi 



(12) 



(13) 



The the sum of the virtual and counterterm corrections, A^^^- and A§^-, is ultraviolet 
finite, as it should be, but still infrared divergent. The calculation has been performed 



in the dimensional reduction scheme [HI, which preserves supersymmetry at the one 



loop level. However, the results are the same in both the dimensional regularization 15 



and dimensional reduction schemes if the quark mass counterterms, which are needed to 
renormalize the ^qq' interaction, are shifted appropriately by a finite amount |jl6[- The 



infrared divergence, which is regulated by introducing a fictitious mass A for the gluon, is 
cancelled after including the real corrections Af^^. 



3.1 Virtual Corrections 

The 0{as) virtual corrections in eq. ([I3|) can be decomposed in the following way, 

A ^ — A^ 4- A™^ 4- A^^ -\- A^ 



(14) 



with A^, A™^, A^^ and A^ the contributions from the diagrams with gluon exchange, 
squark mixing, the quartic squark interaction and the gluino exchange respectively. 
The contribution from the gluon exchange in the $gg' vertex is given by 



A^ 



G 



<S>ij 



Bo (m| , A, mgj + Bq (m?, , A, m^p - 5o (M| , m^-^ , m^^ ; 
+2 (m| + m?, - Ml) Co(m? , M|, m?, ; A, mq^,mq') 



(15) 



with Bq and Cq the Passarino-Veltman scalar two- and three point functions], the ex- 



pressions of which can be found in Ref. [^ 



The contribution from the off-diagonal self energies of the external scalar quarks eval- 
uates to 



A mix 



G 



^ij' 



+ 



m\ - m?, 
G^i'j 



C2e,S2e, (v4o(ml ) - Ao(m| )) + 4:C2e,mg>mgBo{ml,mg, rrig,) 



"q'l 



^h - ^l 



C2e,S2e,{Ao{mlJ - v4o(m?J) + 4c2e,mgm^5o(m|,m^,mg) (16) 



where j' = 3 — j and i' = 3 — i and Aq the scalar one-point function ^7[ . 



The contribution from the loop diagram involving the quartic squark interaction is 
given by 

Ati = - E SiG^uSf^Bo{Mlm^^,m^.) (17) 

fc,/=l,2 

with the matrix S^ defined as 

S^ = ( ^^^' ~'''^ ] . (18) 

Finally, the gluino exchange contribution in the $gg' vertex reads 

Kij = {vsimgivlv^-, + a\a\) + mg{y\v\ - a\a\)) 

-a,{mq{v\a\ + a\v\:) + mg{v\a\, - a^^))} 5o(m|, m^, m^) 
+ [v,[m^[v\v\, + 44,) + mgivy^^ - 44,)) 

+as{mg,{vlai, + a'-vl,) - mg{vlai, - a^-t;^-,))} 5o(m?, ,m^,mgO 
+ {vs{mg{vlvl, + a\a\,) + mq,{v\vl, + a\a\,)) 

-as{mg{v\a\, + a^vl,) - mg,{v\a'-, + a^-t;|,))} ^o(^|, mq, mq>) 
+ [as m-g {v\a\, - a\v\,){Ml - (m^ - mq,f) 
-Vs m-g {v\vl, - a\a\,){Ml - {rUg + mg,f) 
+0's{v\a'~, + a\vl,){ml,mq - m\mq' - (m^^ - mqmq>){mq - niq^)) 

-t^siv'qVl' + a'gai,){ml,mq + m^rriq, - {rrig'^ + mqmqi){mq + m^/))} 

xCo{ml^,Ml,m^^>,mg,mq,mq>) (19) 

Here, Vg and a^ denote the scalar and pseudoscalar couplings of the Higgs boson $ to the 
quarks, which for H, A and H^ are given by 

Hqq : Vg = 2\plmqr\ , a^ = 
Aqq : a^ = 2y2mg(tan/5)"^-^3 ^ f ^ = 
i7+t6 : Vgjas = 2{mb tan/? ± mt ctg/3) (20) 

and f ^, a'g are the reduced gluino-quark-squark couplings which read 



-(J' "-g 



I, . , , o 1 



J = «g- = 2(^^9 ~ *«J ^^^ 4 = -^1 = 2^^s, + sej . (21) 



We have verified that these expressions agree with those given in Ref. |^ in the case 
of the charged Higgs boson decays into tb final states. 
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3.2 Counterterm corrections 



The counterterm corrections consist of the external squark wave-function renormahzation 
and the renormalization of the Higgs-squark interaction. As discussed previously, this 
interaction involves the quark masses, the squark mixing angles and the trilinear squark 
couplings, and a counterterm for each of these parameters will be needed. However, the 
trilinear coupling can be expressed in terms of the quark and squark masses and the 
squark mixing angle 



Aq = — se^ce^ (m\ - m^ - /i(tan/5)" 



-211 



f22) 



with the parameters /i and /3 not renormalized by strong interactions. Therefore, the coun- 
terterm for the trilinear coupling 6Aq is fixed by the quark and squark mass counterterms 
(5mg and (5mg., and by the mixing angle counterterm bQq 



<5A 



2 2 

m|^ — m 



2;e^(^'-.*«' 



S20„ 



bVTir. 



m„ 



+ S20, m^ 



6m, 



91 



5m, 



91 



m„ 



S2e, rrig^ 



12 



■m„ 



(23) 



We perform the renormalization programme in the on-shell scheme where the quark 
and squark masses are defined as the poles of their respective propagators Srrig^ = S^j(m?) 
and SrUq = S?j(mg). The squark wave-function renormalization constants 6Zg- are defined 
as usual in such a way that the residues at the poles are equal to one. 



Finally, as done in Ref . [|T^ , the mixing angle counterterms are fixed by the requirement 
that the renormalized self energy for one of the squarks should remain diagonal, and we 
will choose this squark to be the lightest one. This is similar to the treatment made in 



Ref. 19 



The various counterterms then read 



1 SZ, . (,nl . ,„/ - ,4) B„'(,4, ,„„ ™,) - B„(™|. ™„ ™,) . B„(,4, A,™,) 



+2{-lYs2e, m-gTriq Bo\m\,mq,mg) + 2m| Bo\m\, X,mg^ 



(24) 



mjm,^ = -(m^ + m/-m|)Eo(m|,m„m,)-2m|5o(m|,A,m,J 
-Aoinig) - Ao{mg) + 2 [(^ + 4e,) ^o("^gJ + sj^^ Ao{mg^,) 



-2(-l)*S20, m-gmg Bo(m? ,mg,mg 



(25) 



6ma 



m„ 



2Bi{ml,mg,\) + 4So(mq,mg, A) + Bi{ml,mg,mgJ + Bi{m^g,mg,mg^) 



+S2e, 



rria 



rrir, 



Bo{mg, rrig, m^J - Bo{mg, nig, nig^] 



(26) 



and 



5e„ 



ml - m| 



4mgmgC2e, 5o(m?^, mg,mg) + C2e, S20, (Ao(mgJ - Ao{mq^' 



where again i' = i — 3 in eq. ( P^D and the functions Bq and i?i are defined as 

d 



B'^iq , mi, 1112) = —Bo{q,mi,m2 
B,{q\m,,m2) = ^ 



A{mi) - A{m2) + {m\ -m\- q'^)Bo{q'^, mi, -022) 



The sum of all counterterm corrections is then given by 



ACT 



Ig^JsZ,^ + SZr) + ^ 5m, + ^ Sm, 



dmt 



dm,b 



^ ^*^^- SA, + ^^ 6 A, + ^^ 69, + 



dG. 



mj 



dAt dA, 36, 36, 

for the charged Higgs boson, whereas for neutral Higgs decays it reads 

^%lj = ^G<t,ij(6Zq^ + 6Z, 



'13 



3m,n 



'^ ' 3A, ^'"^ 



36^''^ 



(27) 



(2J 



(29) 



(30) 



Note that only the wave-function counterterm involves infrared divergences and the 
gluon mass A can be set to zero everywhere except in the last term of eq. (|2^) . 

In Ref. 0, the renormalization of the quark and squark masses has been performed in 
the DR scheme, i.e. by subtracting the poles and the related constants in the dimensional 
reduction scheme; in addition there was no explicit counterterm for the trilinear squark 
coupling Aq. The renormalization scheme adopted here is more convenient, since the 
conditions we used for the masses allow for a straightforward physical interpretation. 



3.3 Real corrections 

Finally, the real corrections with an additional gluon in the final state need also to be 
included. In agreement with Ref. p|, we obtain 



aR 



8M^Gij,ij 



AV2(M|,m|,m2,) 



{Ml - ml - m|,)/i2 - mpn - m?,/22 - h - h (31) 



where the phase space integrals lij and lij are defined as 

L i Qj fCq d rCj^ d rCj ['^^ i^i i^j '^Q J 



■■IJ 



'«J 



7r^ 



2Eg 2Ei 2Ej {2kg.h){2kg.kj 



i. Id tXtQ d rCj^ d rhj ['^^ i^i i^j i^n 



TX^ 



2Eg 2Ei 2Ej 



[Ikg.kijj 



(32) 



with kg, kij the four- momenta of the gluon and the squarks qij. In terms of Dilogarithms, 
this gives 



^V_T<Jlio 



\Mq,mq-mq 



^% = ^{(M| - ml - m|) [ - 2 log (^ \^ '^ j log/3o + 2 log^ Po 



- log2 A - log2 /?2 + 2Li2(l - P^o) - Li2(l - P!) - Li2(l - Pi) 
+2Af log J^^^l^^J + 4Af + (2M| + rn| + m|) log/5o 



+ (M| + 2m^p log/52 + (M^ + 2mJ-J log A} (33) 

with Ajj = A(M|,m|,m?,) and 



/3o 



iW| - m? - ml + Xll' 



2mqmq'^ 



Ml - ml + ml - Xlj' Ml + ml - ml - >^ 
A = ^7 , and /?2 = ^^7 (34) 



^j 



4. Numerical Results 

In the numerical analysis of the QCD corrections we have chosen mt = 175 GeV and 
mfe = 5 GeV for the top and bottom masses, s^ = 0.23 for the electroweak mixing angle, 
and as = 0.12 for the strong coupling constant. In Figs. 1-4, ttt,^-^ always denotes the 
mass of the lightest scalar top quark. The Higgs masses are fixed by Ma = 600 GeV, 
and for tan/3 we choose the value tan/3 = 1.6. We will always assume that the left- and 
right-handed scalar masses are equal, m^^ = m^-^ = ^ = ""^6 • 

The dominating decay modes of the neutral Higgs bosons are the decays $ — > tjtj into 
scalar top quarks. In order to visualize the effect of the QCD corrections. Figs. 2 and 
3 display the decay widths of the heavy CP-even Higgs boson H and the pseudoscalar 
Higgs boson A into top squark pairs in Born approximation (solid lines) and including 
the QCD corrections for two values of the gluino mass: m,g = 200 GeV (dashed lines) and 
m.g = 1 TeV (dotted lines). 

Fig. 2a reflects the situation for unmixed top squarks, where the choice At = —fi ctg/? 
yields a diagonal stop mass matrix eq. (^. Fig. 2b refers to maximal mixing 6t ^ —n/A 
in the scalar top sector, with At = 250 GeV. fi = —300 GeV is a common input for 
Figs. 2a,b. The decay widths are significantly larger for the case of mixing, being further 
increased by large QCD corrections up to nearly 50%, whereas in the unmixed case the 
QCD corrections decrease the Born width significantly for the major part of the ti mass 
range. Only close to the phase space boundary, the higher order contribution is positive. 

The non-diagonal decay mode H -^ tit2, shown in Fig. 2c for fi = 100 GeV and 
At = 150 GeV, has a similar signature as the diagonal one in Fig. 2b, but it is suppressed 

10 



by several orders of magnitude and is thus of less interest. 

For the CP-odd A boson, the non-diagonal decay is the only allowed decay mode. The 
width is comparable in size to that of the diagonal H boson decay. We discuss the impact 
of QCD corrections in terms of two examples: Fig. 3a corresponds to the parameters 
H = —300 GeV, At = 300 GeV. Here the corrections show a similar pattern as for the 
diagonal decay mode of the H with maximal mixing, Fig. 2b, having a milder dependence 
on the gluino mass. For other values of the n, At parameters where the lowest order A 
width is smaller, the effect of the QCD corrections can be much more dramatic, as shown 
in Fig. 3b for fi = 100 GeV, At = 250 GeV. Depending crucially on the gluino mass, the 
QCD loop contribution can be either positive {rrig = 200 GeV) or negative {rrig = 1 TeV), 
in both cases up to the order of 50%. The kink in Fig. 3b corresponds to the threshold 
rrit^ = irtg + mj in the ^2 wave function renormalization. 

For completeness, we also present numerical results for the decay of the charged Higgs 
boson H^ -^ tibi in order to demonstrate the impact of the QCD contributions. This 
decay was already studied earlier in the literature using a DR renormalization scheme. 
As for the decay of the heavy CP-even Higgs boson H, we consider the two situations of 
unmixed top squarks (Fig. 4a) obtained by At = —fictgP, A^ = —[i tan/3, and for the 
case of mixing with At = Ai, = -100 GeV (Fig. 4b). Both Figs. 4a,b are for /i = 200 GeV. 
The QCD corrections are large and negative, independent of top squark mixing. They 
can decrease the partial decay width by almost an order of magnitude. 

As seen previously the magnitude of the QCD correction strongly depends on the 
gluino mass. In Fig. 5a we show the dependence of the H boson decay width on nig, 
which is treated here as an independent parameter. The two curves correspond to the 
case of unmixed top squarks (dotted) and stop mixing (dashed-dotted) with At = 250 
GeV. The other parameters are fixed according to // = —300 GeV, m^^ = 150 GeV 
and the solid and dashed lines indicate the respective widths in the born approximation. 
Whereas for the unmixed case the QCD contributions decrease the width for large rrig, 
the decay width increases with rUg for the mixing case. The asymptotic dependence on 
rrig is logarithmic; all terms linear in rrig cancel in the final result. Such a logarithmic 
behavior has been first observed and discussed in the case of the decay of squarks into 
quarks and photinos [^ . 



The asymptotic dependence on rrig can be cast into a simple form. In the notation 

and normalization of eqs. (7-10), the decay amplitude G^tih ^^^ ^^^ process A — > tit2 
gets an additive gluino mass dependent contribution of the type 



rrit 



V2 



^GAhh - -m iA-Mi^)ctgf3-fi -^\og{mg) 



4a 



37r 

rritH 4 0^ 

log{mg) . (35) 



For the H boson decays, the structure is shghtly more complicated, but still very compact: 
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The coupling matrix Ghe ^^ eq. (|^ gets the following additive contribution: 



4 n 

^Ghu = -^log{m-,)ACHu (36) 

6 n 



with the matrix 

^Cjjti = 2{R) ChuR 




.t ^2 



mtHAt + Ml^yi+fir'^] 



mt[{At + Ml^)r\ + fi rl] 8r\mf 



K (37) 



where C^n is the matrix in eq. (|^). For the special case H —>■ titi we obtain: 

~ 2CHi,t, + -^{Sm^ctg/? + 2se,ce, mt[{At + Mi^)ctg/5 - /i]| (38) 

From this expression, the different behaviour of the m^-dependence for the mixed and 
unmixed case in Fig. 5 can be qualitatively understood. It should be noted that a strong 
rrig-dependence for intermediate values of irig is introduced via the finite part of the vertex 
diagram with q and g in the internal lines (Fig. lb). It decouples for large m^, but the 
decouphng is very slow. 

This gluino vertex diagram is also responsible for a finite pseudoscalar decay width in 
case that the tree level coupling of A to gig2 is zero, i.e. for correlating /x and Aq such that 
the entries in eq. (P) vanish. Together with the counterterm from the renormalization of 
the trilinear coupling Aq, eq. (P5|), the 1-loop amplitude is ultraviolet and infrared finite 
and gives rise to a loop-induced finite decay width. This situation is shown in Fig. 5b 
for y4° — »• tit2 for the parameters At = ^ tan/3, /i = ±100 GeV and m^^ = 200 GeV as a 
function of the gluino mass. The strong dependence on Tjig is to a large extent due to the 
finite part of the gluino vertex diagram, eq. (|19|). 



5. Squark decays into Higgs Bosons 

For completeness, let us also discuss the decays of scalar quarks into lighter squarks and 
Higgs bosons. In practice this situation can occur when there is a large mass splitting 
between the heaviest and the lightest top squarks where the decay 

h^h + h/A/H (39) 

can occur, and between the stop and sbottom squarks where the decays 

ts ^ H+ + bi or bi^ H- + h (40) 

can be kinematically allowed. In the previous equation, hi is the lightest bottom squark 
which is identical to the left-handed h squark in case of no mixing in the sbottom sector; 
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decays of the heaviest bottom squark into neutral Higgs bosons can also be possible for 
large tan/5 values where the 61-&2 mass splitting and the sbottom mixing angle can be 



large enough. These decays have been discussed at tree-level in Ref. |T0[ for instance. 



The partial decay width of a squark g^ to a Higgs boson $ and a lighter squark q'j is 
given by a relation similar to eq. (^ but without the color factor 

r°(g. - $5) = — £^ V/^K,M|,m|,) ( G^., f (41) 

m 

with the coupling G$jj, with $ = i7, A and H^ given as in section 2. The couplings of 
the lightest CP-even Higgs boson h to squark pairs can be obtained from those of the H 
boson by simply replacing a by a — 7r/2. 

The QCD corrections to the decay width are given by the Feynman diagrams shown 
in Fig. 1 with now g^ and $ being in the initial and final state respectively. The QCD 
corrected decay width can be written as 

T\q, -. $5) = ^^^Ai/^(m|,M|,m|,) G^,, A^,, (42) 

Q\/2Txm% TT ^' ^J 

Hi 

where, similarly to eq. (0), the correction factor /\^ij is given by 

A^,,- = A^,^. + AgJ. + Af,^. . (43) 

The virtual corrections A^^^ and the corresponding counterterm A^^- are given by the 
same expressions as in section 3. In the real corrections, however, the role of qi and the 
Higgs boson $ have to be interchanged: 



,2 









(M^ - mi - mi ) /02 - mi /qo - mi I22 - h - h (44) 



where the functions Iik |^ have arguments 

Ilk = Iik{mg^,mg>,M^, A) (45) 



^j 



Fig. 6 shows the partial decay widths of the heaviest top squark ^2 into its lighter 
partner ii and the lightest neutral CP-even Higgs boson h (6a) and the pseudoscalar 
Higgs boson A (6b) as a function of the lightest top squark mass. Again, the solid 
lines show the tree-level decay widths, while the dashed and dotted lines show the QCD 
corrected widths for m^ = 200 GeV and rUg = 1 TeV, respectively. To have the t2-ti 
mass splitting large enough while leaving the decay rates sizeable, we relax the condition 
m,f^ = m^i^ for the stop mass parameters. We chose /x = —300 GeV, tan/3 = 1.6, At = 300 
GeV and mr = 500 GeV; mr is then varied from 200 to 430 GeV. 

In Fig. 6a, the pseudoscalar Higgs boson mass is taken to be Ma = 400 GeV, leading 
to a value M^ — 75 GeV. This value slightly varies with m,];^, since we include the large 
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radiative corrections in the Higgs sector which increase with mf and with the logarithm 
of mi^mi^. The ^2 -^ hti decay width is about 0.1 GeV for small ti masses and the QCD 
corrections increase the decay width by an amount of (9(10%). 

In Fig. 6b, the pseudoscalar mass is chosen to be Ma = 100 GeV. The ^2 -^ ii^ decay 
width is an order of magnitude larger than in the previous case, but the QCD corrections 
are now negative, decreasing the partial widths by 20%. Again, the kinks in Figs. 6 are 
due to the opening of the 1 2 — *> t + ^ threshold. Above this value, this decay becomes by 
far the dominant one compared to the decays into Higgs bosons. 

6. Conclusions 

We have calculated the SUSY QCD corrections to the decays of the MSSM heavy neutral 
and charged Higgs bosons into scalar quark pairs, including the mixing in the squark 
sector. A special attention has been paid to the case of stop and sbottom decays which 
can be the dominant decay channels in a large area of the MSSM parameter space. The 
QCD corrections turn out to be quite substantial, enhancing or suppressing the decay 
widths in Born approximation by amounts up to 50% and in some cases more. The QCD 
corrections depend strongly on the gluino mass; however, for large gluino masses, the 
QCD correction is only logarithmically dependent on rrig. Contrary to the case of Higgs 
decays into light quarks, these QCD corrections cannot be absorbed into running squark 
masses since the latter are expected to be of the same order of magnitude as the Higgs 
boson masses. 

We have also calculated the QCD corrections to the decays of heavy squarks into lighter 
squarks and Higgs bosons. Cases of interest are for instance the channels ^2—^1^1 + h/A. 
The QCD corrections are at the level of a few ten percent and can enhance or suppress 
the tree-level decay widths. In both cases, Higgs decays into scalar quarks and squark 
decays into Higgs bosons, the QCD corrections are therefore important and should be 
included when discussing these decays. 

After completion of this work, another paper appeared [^ which deals with the 



decays of the MSSM Higgs bosons into squark pairs. It is an extension of the earlier 
work for the H^ decay, but the discussion is performed also in the on-shell scheme. 
The difference to our renormalization scheme is only a slightly different treatment of 
the squark mixing angle renormalization, which results in a minor finite shift in the 
mixing angle counterterm. We have checked that all analytical expressions agree with 
ours. Numerical comparisons also show good agreement; in our case we show that the 
correction can strongly decrease the decay widths and determine the origin of the large 
corrections. The decays of heavy squarks to lighter squarks and Higgs bosons has not 
been discussed in Ref. |]22[| . 

Acknowledgments: We thank Roland Hopker for his collaboration in the early stage 
of this work. 
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Figure Captions 

Fig. 1: Feynman diagrams relevant for the 0{as) QCD corrections to the decay of a Higgs 
boson $ into scalar quark pairs, (a): gluon exchange, (b) gluino exchange, (c) mixing 
diagrams, (d) quartic squark interaction, (e) self-energy and vertex counterterms, 
(f) real corrections. 

Fig. 2: Partial widths [in GeV] for the decays of the neutral CP-even Higgs boson into 
top squark pairs, H -^ Utj, as a function of of tti^^. The Higgs mass is fixed by 
Ma = 600 GeV and tan/5 = 1.6; ^ = -300 GeV, A = -/ictg/5 (a); /i = -300 GeV 
and At = 250 GeV (b); /i = 100 GeV and At = 150 GeV (c). The solid lines are 
for the partial widths in the Born approximation, while the dashed and dotted lines 
are for the partial widths including QCD corrections for rrig = 200 GeV and 1 TeV 
respectively. 

Fig. 3: Partial widths [in GeV] for the decays of the neutral CP-odd Higgs boson into top 
squarks, A -^ iii2, as a function of m^-^. The Higgs masses are fixed by M^ = 600 
GeV, tan/5 = 1.6; /i = -300 GeV and A = 300 GeV (a); /x = 100 GeV and A = 250 
GeV (b). The solid lines are for the partial widths in the Born approximation, while 
the dashed and dotted lines are for the partial widths including QCD corrections 
for rrig = 200 GeV and 1 TeV respectively. 

Fig. 4: Partial widths [in GeV] for the decays of the charged Higgs boson into the lightest 
top and bottom squarks H^ -^ tibi as a function of m^-^ . The Higgs mass is fixed by 
Ma = 600 GeV and tan/5 = 1.6; /i = 200 GeV, A = -/ictg/3 and A = -/x tan/3 
(a) and fi = 200 GeV and At = Af, = —100 GeV (b). The solid lines are for the 
partial widths in the Born approximation, while the dashed and dotted lines are 
for the partial widths including QCD corrections for rrig = 200 GeV and 1 TeV 
respectively. 

Fig. 5: Partial width [in GeV] for the decay of the neutral CP-even Higgs boson H into 
the lightest top squarks (a) and of the neutral CP-odd Higgs boson A into top 
squarks, as a function of the gluino mass. The Higgs boson masses are fixed by 
Ma = 600 GeV and tan/5 = 1.6. In (a) we take /i = -300 GeV, m^^ = 150 GeV 
and At = -/ictg/3 (dotted) and fj, = -300 GeV, rrif^ = 150 GeV and At = 250 GeV 
(dashed-dotted). The solid and dashed lines corrrespond to the respective born 
widths. In (b) : At is fixed to make the lowest order width vanish. At = fi tan /5 and 
the stop masses are fixed by m^-^ = 200 GeV; n = 100 GeV (solid) and fj, = —100 
GeV (dashed). 

Fig. 6: Partial decay widths of the heaviest top squark t2 into ti and the lightest neutral 
CP-even Higgs boson h (a) and the pseudoscalar Higgs boson A (b) as a function of 
rrit-^ . The solid lines show the tree-level decay widths, while the dashed and dotted 
lines show the QCD corrected widths for rrig = 200 GeV and 1 TeV, respectively. 
H = -300 GeV, tan/5 = 1.6, A = 300 GeV and m^^ = 500 GeV while m^-^ is 
varied from 200 to 430 GeV. For the Higgs masses we take Ma = 400 GeV (a) and 
Ma = 100 GeV (b). 
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